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Abstract 

We investigate the implications of non-standard interactions on neutrino oscillations in the 
OPERA experiment. In particular, we study the non-standard interaction parameter e^r- We 
show that the OPERA experiment has a unique opportunity to reduce the allowed region for this 
parameter compared with other experiments such as the MINOS experiment, mostly due to the 
higher neutrino energies in the CNGS beam compared to the NuMI beam. We find that OPERA 
is mainly sensitive to a combination of standard and non-standard parameters and that a resulting 
anti-resonance effect could suppress the expected number of events. Furthermore, we show that 
running OPERA for five years each with neutrinos and anti-neutrinos would help in resolving the 
degeneracy between the standard parameters and e^r- This scenario is significantly better than 
the scenario with a simple doubling of the statistics by running with neutrinos for ten years. 
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I. INTRODUCTION 



Neutrino oscillation physics has definitively entered the era of precision measurements of 
the fundamental neutrino parameters such as the neutrino mass squared differences {i.e., 
Amg^ and Am^i) and the leptonic mixing parameters {i.e., 612, ^13, 623, and 6). In partic- 
ular, the Super-Kamiokande, SNO, KamLANDjK2K, and MINOS experiments have given 
valuable information on these parameters Q, [3, 3, 0, [s, 6]. 

The precision measurements open up the possibility to investigate if neutrino flavor tran- 
sitions are governed by neutrino oscillations only or if they are, in the next-to-leading order, 
a combination of neutrino oscillations and some other new physics mechanism. However, to 
leading order, there exist clear evidences that neutrino oscillations constitute the underlying 
physical model for neutrino flavor transitions. The next-to-leading order mechanism could 
e.g. be non-standard interactions (NSIs), mass varying neutrinos, neutrino decay, neutrino 
decoherence, etc. or some combination thereof. 

In this work, we will study NSI effects at the OPERA experiment Q], which is an ex- 
periment that consists of a massive lead/emulsion target (the OPERA detector) located at 
LNGS in Gran Sasso, Italy, receiving its neutrino beam, originally consisting almost ex- 
clusively of Ufj^, from CERN in Geneva, Switzerland. The baseline length is approximately 
732 km and the CNGS z/^ beam has an average neutrino energy of E^, ^ 17 GeV. The 
OPERA experiment is especially designed to observe events from the z/^ neutrino 
oscillation channel. In fact, no previous experiment has investigated this channel or ob- 
served neutrinos of a different flavor than that originally produced at the source (although 
the neutral-current measurements at SNO imply that solar have oscillated into a different 
flavor). Thus, the OPERA experiment presents a unique opportunity to study direct ap- 
pearance of Sj. In this work, we will not try to describe the origin of the NSIs, but adopt 
a purely phenomenological point of view. In particular, NSIs can modify the production, 
the propagation in matter as well as the detection of the neutrinos. We will concentrate on 
the simplifled scenario in which NSIs only affect the neutrino propagation. 

Previously, investigations of NSIs that are of importance for this work have been presented 
in the following papers: In Ref . [sl , a two-flavor neutrino analysis of the so-called atmospheric 
neutrino anomaly has been performed, which effectively bounds the NSI parameters in the 
H-T sector, e^r — ^ and Err — ^' ■, to —0.03 <e< 0.02 and \e'\ < 0.05 at 99.73 % confldence 
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level. Although these bounds may seem quite restrictive, it has been shown that at least 
the bound on e-rr is severely weakened when considering the full three-flavor framework 
(allowing Err to be of 0{1) or larger, depending on the values of See and Ser As will 

be shown later in this work, the limit that could be put by the OPERA experiment would 
be insensitive to whether the two- or three-flavor scenario is studied simply because of the 
relatively short baseline. In addition, in Ref. the authors have come to the conclusion 
that it would be possible to observe NSI effects at the OPERA experiment (and the ICARUS 
experiment) if e^,- > (9(10~^). Next, in Ref. 12], the Kamioka-Korea two detector setup 
has been investigated, which could also give restrictions on the NSI parameters and 
Err- Recently, in Ref. 13|, a study of the OPERA experiment (in combination with the 
MINOS experiment) has been presented with the conclusion that it is not very sensitive to 
the NSI parameters Eer and Err- However, it was found that the v-r sample is too small to 
be statistically significant to improve the limits on the NSI parameter Err- Nevertheless, 
this analysis did not include a study of the relevant e^,- which, due to the energies and the 
baseline involved in the OPERA experiment, is the only NSI parameter appearing to leading 
order in L in the Va — ^ i't flavor transition. 



In general, neutrino oscillations and NSIs in terrestrial neutrino experiments have 
studied extensively in the literature, using the neutrino factory project 



been 
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2l| and other different neutrino facilities (like super-beams and /3-beams) [2 
261] to assess the impact of the NSI effects in neutrino physics. 
This work is organized as follows. In Sec. [Tll we will present analytic considerations for 
the NSIs that we assume for the OPERA experiment. In addition, we will comment on a sort 
of anti-resonance effect that is in the vicinity of being detectable in the OPERA experiment. 



Next, in Sec. Illlt we will give the numerical setup with the GLoBES software 27|, |28| that 
we use for our simulations of the OPERA experiment. Then, in Sec. llVt we will show our 
numerical results for the OPERA experiment using GLoBES. Finally, in Sec. |Vl we will 
present a summary of the work as well as our conclusions. 
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II. ANALYTIC CONSIDERATIONS 



We consider effective non-standard interactions of the form 



^NSi = -^ E ^S[/7ni + «7')/][^7.(l-7')^/.], (1) 

f — u,d,e 

a=±l 

where / is summed over the matter constituents and the parameters e:^^, which are the 
entries of a Hermitian matrix e-^"", give the strength of the NSIs. In a manner completely 
analogous to the derivation of the normal matter effect, these interactions will result in an 
effective addition 



-f^NSI = V 



^ ^ee ^efj, ^er 



^efi ^fJ-fJ. ^/JT 



(2) 



fa 



to the neutrino oscillation Hamiltonian in flavor basis, where Ea/s = Ylf a^ap-^f /-^e and 
V = ^/2GfNi,. Notice that, apart from the bounds on e^r and Ett given in the Sec. [H we are 
not aware of any paper discussing direct bounds on the effective par ameters Ea/s- However, 
experimental limits on the parameters can be found in Refs. 29|, [30 1, which imply that 
ke/Jl — C^(10~^) and \eI^\ < 0{1) [20]. Thus, we can assume that the effective parameters 
Eee and Ee^ are bounded at the same order of magnitude as their corresponding parameters 

The full three-flavor Hamiltonian describing neutrino propagation in matter is given by 

H = diag(0, Am^i, Aml^)U^ + Husw + H^si, (3) 

where U is the leptonic mixing matrix, Amfj = mf — rrip and Husw is the addition from 
the standard matter effect. Due to the quite large neutrino energy E^, = (9(10) GeV and the 
relatively short baseline L ~ 732 km, both AmliL / {2E^) <^ 1 and VL <C 1, where V is the 
matter potential ~ 1.1- 10"^'^ eV in the Earth's crust (p 2.7g/cm^) fsi)]. Thus, neutrino 
oscillations will not have time to fully develop. As a consequence, the main characteristics 
of the flavor transition probabilities will be given by truncating the flavor evolution matrix 
S = exp(— iifL) at order L, resulting in 

^ ~ 1 - iHL. (4) 
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The off-diagonal neutrino transition probabilities are then given by 

Pa,3 = l^/j^r ^ \H^aL\^. (5) 

The diagonal neutrino survival probabilities in this expansion are given by the unitarity 
condition P^a = 1 — S/j^a ^ajS- As can be observed from this consideration, the transition 
probabilities will only be affected by the corresponding NSI element {i.e., P^p just depends 
on the NSI element Spa)-, while the survival probabilities depend on the two off-diagonal NSI 
elements associated with the flavor {e.g., P^^ is affected by Sefj. and e^r)- As expected, the 
diagonal NSI parameters do not enter at short baselines. Clearly, this is not true in general 
and at higher orders in L, where the NSI parameters will enter all of the neutrino oscillation 
probabilities. As an example, the NSI parameter e^t will enter the flavor evolution matrix S^r 
at 0{L'^) and then to 0{L^) in the transition probability P^j. (unless there is no interference 
between the L and terms). From the above consideration, we can conclude that the NSI 
parameter of most interest for the OPERA experiment is e^,-- That the parameters Eer and 
Err are not important has been already shown in Ref. [3]. 

The main physics goal of the OPERA experiment is to actually observe oscillations of z/^ 
into Ur- With the effects of E^r included, the transition probability P^,- is given by 



P — I ^ I 

/XT I '-'TfJ. I 



4 sin(2M^ + <.^' 



L' + 0{L'), (6) 



where we have neglected the small mass squared difference Am^i. From this consideration 
follows that there is a degeneracy between the standard neutrino oscillation parameters and 
the NSI parameter e^^ as scenarios with the same value of jc^g sin(26'23)Am3^/ {4:E^) + e*^^V\ 
will lead to the same neutrino oscillation probability. Even if the degeneracy is broken 
by the energy dependence of the first term, we still expect some parameter correlations 
when analyzing the outcome of an experiment. It is also interesting to note that the 0{L'^) 



contribution to P^,- vanishes when 



<. = -4^sin(2M (7) 



simply due to the fact that Srfi = in this case. The condition clearly shows that this 
can happen only for real e^^- We will use the term anti-resonance to refer to this scenario 
as it, in some sense, is the opposite of the MS W- resonance: in the standard picture of 
neutrino oscillations, the matter effects cancel the difference between the diagonal terms 



and the effective mixing angle is maximal, whereas in the situation with NSIs, the matter 
effects cancel the off-diagonal terms and the effective mixing angle is minimal {i.e., zero). 
In a pure two-ffavor scenario, the anti-resonance is valid to all orders, while transitions can 
be induced to higher order in L by other off-diagonal elements in the case of three-ffavor 
oscillations. For the peak energy of -Ej, ~ 17 GeV in the CNGS beam, the anti-resonance 
would occur for e^,- ~ —0.3 with the result that no Ur events would be observed. Note that 
a similar conclusion applies in the case of inverted mass hierarchy, from which e^r — +0.3 
if Amg^ — > — Arrig^ (neglecting the small effect of Am2i). This also applies to the case of 
anti-neutrinos, where we have V —V and Sais — > e^^^. In both cases, this also gives an 
estimate of the order of magnitude of the NSIs that OPERA will be sensitive to, as the 
expected number of z/^ events is low. 

Finally, we want to mention that a similar effect could exist in the — »• z/g transition. 
In fact, 

2 

L2 + 0(L3), (8) 

where 6 is the standard CP-violating phase in the unitary leptonic mixing matrix, a = 
Am2i/Am|^ is the ratio between the mass squared differences, and we have neglected a 
term proportional to Si^a. In this case, the external bounds on Sefi are so stringent that the 
term proportional to a is known to be larger. Thus, an anti-resonance in this channel could 
only be due to an interplay between the two standard terms if 5 = tt. 



lie 



\s, 



e/i I 



A TOO ^ 

[s23 sm{2eis)e'' + acssCis sin(20i2)^ 
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III. NUMERICAL SETUP 



The numerical simulations of the OPERA experiment were performed using the GLoBES 



software 



27 



281], which was extended in order to accommodate the inclusion of NSIs through 
the Hamiltonian presented in Eq. ([2]) with e^e = £efi = Ser = 0. The neutrino propagation 



in matter was then described using the full three- flavor Hamiltonian in Eq. In addition, 
the Abstract Experiment Definition Language (AEDL) file, used to describe the OPERA 



33l | . Unless stated otherwise. 



experiment, was based on the results presented in Refs. [tI. H. 
we have assumed a running time of five years with 4.5 ■ 10^^ protons on target per year, 
in accordance with the OPERA experimental setup, and an effective mass of 1.65 kton [7| . 
Furthermore, the neutral- and charged-current cross-sections were taken from Refs. [3, Q, 



012 = 34.4° ± 1.7° 


Am|i = (7.59 ± 0.21) • 10"^ eV^ 


6*13 = 4.8° ± 2.9° 


Am|i = (2.4 ± 0.15) • lO'^ eV^ 


023 = 45° ± 3.8° 


(5 = 7r/2 



TABLE I: T/ie simulated values of the standard neutrino oscillation parameters and the correspond- 
ing la priors used in the simulations. The central values of the parameters 0\2 and Am2i were 
inspired by the results of the KamLAND experiment [6], whereas the central values of the other 
parameters were inspired by Ref. [37]. We fixed the value of the CV -violating phase to 7r/2, with 
no consequences on our results for the u^- channel. 



35|. The CNGS neutrino spectra are substantially different from zero in the interval between 
1 GeV and 30 GeV (with a peak around ^ 17 GeV). Thus, we divided the signals and 
the corresponding backgrounds into 29 equally spaced energy bins, having checked that the 
numerical results are stable if the number of energy bins is above the order of 10. For the 
baseline length of the CNGS setup (approximately 732 km), the matter density profile was 
assumed to be constant and equal to the value at the Earth's crust, i.e., p = 2.72 g/cm^ 
(or V = 1/1900 km~^) 36|. In all simulations, we have used a full three-flavor neutrino 
framework with central values and la errors of the standard neutrino oscillation parameters 
as given in Tab. [B Normal mass hierarchy, i.e., Am|^ > 0, has been assumed if not stated 
otherwise. 

Regarding the NSI parameters, we performed numerical simulations with different sim- 
ulated values, also taking into account the effects of possible CP-violating phases of the 
non-diagonal entries of the Hamiltonian in Eq. ([2]). The priors set on the NSI parameters 



are chosen according to Ref. 29|], except from Err, which has further been constrained using 
atmospheric neutrino data [10]. As a comparison, we also included the MINOS experiment, 
able to probe the z/^ transition channel, in our simulations. As already mentioned 

in Ref. |lj|, different L/Ey could in general be very useful in order to further constrain 
some of the parameters of Eq. ([2]), since the relative importance of the standard and non- 
standard parts of the Hamiltonian is energy dependent. Our numerical setup of the MINOS 
experiment follows that used in Ref. 24l |. 
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FIG. 1: Sensitivity for e^r <it 95 % confidence level (2 d.o.f.) of the OPERA and MINOS experi- 
ments in the case of no NSIs (the input values of the various e^p = 0^. 

IV. NUMERICAL RESULTS 

In this section, we present the numerical results on the physics reach of the OPERA 
experiment in constraining the new physics parameters Eap- In all figures, we have combined 
both the and — > channels for the OPERA experiment. 

The results have been obtained by marginalizing over the parameters Amg^^ and (if 
not stated otherwise), while keeping the parameters Am^i and 612 fixed, since they are 
irrelevant for the i/^ transition in the OPERA experiment. In addition, the parameter 

6^13 was fixed, since it does not affect the results. We also observed that eee,£eij,,£fifi, and ^er 
do not affect the results, which means that they are fixed to zero in the rest of the work. 

First, in Fig. [T], we present the sensitivity reach for with the OPERA experiment in 



combination with the MINOS experiment 



for a discussion on the sensitivity reach for 



and Err for the same combination, see Ref. 13(]). As can be observed in this figure, OPERA 



is far more sensitive to e^r due to the higher neutrino energy than that in MINOS, which 
can therefore only marginally improve the sensitivity. Thus, in the following we will only 
consider the bounds which can be placed from OPERA itself. 
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FIG. 2: T/ie left panel shows the NSI sensitivity of the OPERA experiment alone for marginalized 
standard neutrino oscillation parameters. The confidence levels (2 d.o.f.) are 90 %, 95 %, and 
99 %, respectively. The sensitivity contours for fixed standard neutrino oscillation parameters only 
differ slightly from this result. The right panel shows how the sensitivity to e^r changes depending 
on the Err prior. The sensitivity contours in the right panel are at 95 % confidence level (2 d.o.f.) 
and the Err priors are at la level. 

Note that, as expected, a similar situation is also valid when considering the detection 
of the z/^ Ue transition in OPERA. In fact, as already stressed in Sec. [Tll Sfir appears to 
leading order in P^^. [Eq. ([6])], but it is subleading in P^g [Eq. ([8])]. We verified that the 



inclusion of the z/„ Up channel does not affect the results on e 



A. Marginalization of Ett 

Figure [2] shows the OPERA sensitivity in the le^rl-f^rr plane (left panel) as well as the 
impact on the e^^r sensitivity given different priors on Err (right panel). As expected, the 
impact of Err is small as long as any reasonable prior is put. This is naturally related to the 
fact that Err enters only in higher order in the oscillation probability P^^. From the left panel 
of the figure, it is evident that Err has to be of (9(10) to significantly alter the prediction 
for E^r- In the rest of this work (including Fig. [1]), the prior put on Err is l^rrl < 1-9 (la 
confidence level). 
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FIG. 3: The sensitivity of OPERA to the NSI parameter e^r d.o.f.) for five (colored regions) 
and ten (curves) years of running time. The simulated value of e^r is zero and the confidence 
levels are 90 %, 95 %, and 99 %, respectively. The diamond corresponds to the simulated value of 



e^r = 0. 



B. Sensitivity to s 



/XT 



In Fig. [31 we show the predicted sensitivity of OPERA to the NSI parameter e^r- It 
is clear from this figure that the sensitivity contours extend in the direction where the 
number of Vr events is constant (basically a circle centered at e^r — —0.3). Furthermore, 
the sensitivity change if running the experiment for a longer time has been indicated. As 
can be seen, this would slightly improve the projected sensitivity. However, running the 
experiment with reversed polarity could significantly improve the sensitivity, see Fig. HI 
Due to the change oiV ^ —V and 6^13 — > s*^p when considering anti-neutrinos instead of 
neutrinos, it follows from Eq. ([6]) that the sensitivity contours for the reversed polarity will 
extend in a different direction than those of the original polarity as can be observed in the 
figure. Thus, the combination of the two polarities could aid in resolving the degeneracy. 
That the anti-neutrino run by itself produces slightly larger sensitivity contours is mainly 
due to the lower cross- sect ion. Furthermore, the figure shows the effects of having different 
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FIG. 4: The sensitivity of OPERA for five years of running time in neutrinos (dark curves) and 
anti-neutrinos (light curves) as well as the combination thereof (colored regions). The sensitivity 
levels correspond to confidence levels (2 d.o.f.) of 90 %, 95 %, and 99 %, respectively. 

simulated values for e^^. 

In Fig. [5l the effects of fitting the data to the wrong neutrino mass hierarchy are shown. 
Again, we can observe that the sensitivity contours extend in the direction of a constant 
number of events, i.e., the circle centered at e^r — —0.3. In the case of e^r = —0.3, the 
circle has radius zero and the allowed region is relatively small. This corresponds to the anti- 
resonance case, where no events are expected. The results of the fit using the wrong neutrino 
mass hierarchy is a simple mirroring of the result with the correct hierarchy, —Sf^r- 
This can be easily understood from Eq. ([6]), where a sign change in Sf^r exactly cancels 
the sign change in Amg^^ associated with changing the neutrino mass hierarchy. Strictly 
speaking, the sensitivity of OPERA to e^r is the union of the sensitivities obtained when 
fitting each mass hierarchy separately. 
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FIG. 5: The predicted sensitivity contours of OPERA for different simulated values ofe^r (2 d.o.f.). 



The simulated values chosen are e^r = (upper-left panel), e 



-0.3 (upper-right panel) , e^r 



—0.6 (lower-left panel), and e^r = —0.3(1 -|- i) (lower-right panel). The colored regions correspond 
to a fit using the correct neutrino mass hierarchy, while the curves are the regions obtained with 
a fit using the wrong neutrino mass hierarchy. The confidence levels are 90 %, 95 %, and 99 %, 
respectively. This figure only includes the results of running for five years with neutrinos. 

V. SUMMARY AND CONCLUSIONS 



We have studied NSIs in connection with the OPERA experiment. Unhke in the previous 
work by Esteban-Pretel et al. ISj, where the focus was on the effective NSI parameters Ser 
and due to the external bounds on the other parameters, we have focused on the NSI 
parameter e^^. The reason for this is that e^r is more important for the — > oscillation 
probability in OPERA due to the relatively short baseline, as can be seen in our analytic 
considerations. 

We have found that OPERA is actually sensitive to a combination of standard and 
non-standard parameters, which can be easily observed in Eq. ([6]). The degeneracy in the 
parameter space, where this combination is constant, is somewhat broken by the energy 
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dependence of the standard term. A much better determination of e^r can be obtained if 
we consider a 5+5 year neutrino-antineutrino run, especially if compared with a 10 year of 
data taking with neutrinos only, the main reason being the different correlations between 
standard and non-standard parameters, visible from Eq. (jH]) with the replacement V —>■ —V 
and e^r ^ ^l^r- 

We have also observed that the uncertainty on e^r can be strongly worsened due to our 
ignorance in the sign of the large mass squared difference Am^^, resulting in a reflection of 
the allowed region for e^,- in the imaginary axis. This holds true in the case of running in 
neutrinos only as well as running with both neutrinos and anti-neutrinos. 

By means of the simple result of Eq. we found that an anti-resonance occurs when 
the standard and NSI parameters cancel. In this case, the effective Hamiltonian element 
H-rfj, vanishes and no z/^ — > Ur events would be observed. 

The above analytic considerations were illustrated by our numeric simulations using a 
modified version of the GLoBES software. In particular, it is evident from Fig. H] that 
running the OPERA experiment for five years in each polarity would be much more efficient 
in constraining the parameter space than running for ten years with neutrinos only. 

Finally, we again want to mention that atmospheric neutrino experiments put constraints 
on e^T- which are better than what OPERA is sensitive to However, these constraints 
have been computed in a pure two- flavor framework without the interference of Eer- It is 
known that the bounds on Err from similar considerations are significantly weakened when 
extending to a full three-flavor framework. Thus, OPERA will provide a very clean and 
complementary bound, since only Efj,r enters into the leading term of the neutrino oscillation 
probability P^r due to the short baseline. 
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